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Abstract. We propose a hierarchy of nonlinearly dispersive generalized Korteweg-de Vries (KdV) evolution equations based on a 
modification of the Lagrangian density whose induced action functional the KdV equation extremizes. It is shown that two recent 
nonlinear evolution equations describing wave propagation in certain generalized continua with an inherent material length scale are 
members of the proposed hierarchy. Like KdV, the equations from the proposed hierarchy possess Hamiltonian structure. Unlike 
KdV, however, the solutions to these equations can be compact (i.e., they vanish outside of some open interval) and, in addition, 
peaked. Implicit solutions for these peaked, compact traveling waves (“peakompactons”) are presented. 
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1. Introduction 

The Korteweg-de Vries (KdV) equation, originally derived to describe long waves over shallow water 
mm, has emerged as a paradigm in the theory of solitons and integrable systems following the pioneering 
work of Zabusky and Kruskal Q. Beyond water waves, KdV arises in numerous other physical contexts 
||4l|5l. Additionally, in modeling various physical phenomena beyond long waves over shallow water, many 
generalizations and modifications of the KdV equation have been derived over the years, including but not 
limited to: the modified KdV equation |0|4l, the KdV-Burgers equation Q, the cylindrical KdV equation 
iTTl . the extended KdV equation ||8l, the Camassa-Holm equation Q, theK{n,m) equations ifTOll . theK*{l,p) 
equations im, the KdV-Kuramoto-Sivashinsky-Velarde equation ifT^ . the generalized integrable KdV 
equation ifT^ eq. (4)], KdV-based higher-order equations ifldlfTSl . hierarchical KdV equations ifT^fTTlfT^ . 
the Destrade-Saccomandi equation ifT^ eq. (17)], the Jordan-Saccomandi equation ll20l eq. (4.14)] and 
K{n,m) equations with non-convex advection 1211. 

In the present work, we propose another hierarchy of generalized KdV equations based on modifying 
the Lagrangian density whose induced action functional is extremized by the KdV equation. Furthermore, 
we show that several members of this proposed hierarchy of nonlinearly dispersive evolution equations 
have been recently derived to describe wave propagation in generalized continua (of the type proposed by 
Rubin et al. (221) with an inherent material length scale. Finally, peaked, compact traveling wave solutions 
(“peakompactons”) of equations in the proposed hierarchy are constructed. 


2. A hierarchy of nonlinearly dispersive generalized KdV equations 


Let (p = (p{x,l) be a field. In the present work, motivated by the form of the Lagrangian density whose 
corresponding action functional is extremized by the KdV equation (see, e.g., f23[ eq. (2)]), we consider the 
class of Lagrangian densities given by 
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where subscripts denote differentiation with respect to an independent variable. For convenience, we have 
taken n,m>0 and, below, n and m typically take on integer values. However, the derivations below can be 
easily modified fo hold for any real n and m. The corresponding Hamiltonian density can be obtained by the 
Legendre transformation (see, e.g., (241 §18]): 
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The corresponding action functional / AxAt ^ is extremized by requiring that its first variation vanish 
(see, e.g., (24l §11 and §35] for the details), namely 

d d 

d(p dt \ dq)t J dx \ dq)xj dx^ 

Substituting eq. ([B into eq. ([3]l, assuming (p is smooth enough so that cpxt = (ptx and rearranging, we obtain 
the partial differential equation (PDE): 

(pxt + (<p.v)>xv + [(<P«)'"]„ = 0- (4) 

Interpreting (p as a potential, let u{x,t) = (px(x,t), then eq. dBl becomes 

Mf + u”Ux + {{Ux)"']xx — 0) (5) 

giving the canonical form of the proposed hierarchy of generalized KdV equations as parametrized by {n,m). 



2.1. Conserved quantities 


One standard approach to constructing conserved quantities is to restrict to localized waveforms, specif¬ 
ically such that (p ^ (p± “sufficiently fast” as x —)• ±oo, where (p± are constants. In turn, this condition 
implies that m —>■ 0 “sufficiently fast” as x —> ±oo. Then, multiplying eq. ® by for some integer k, 
integrating overx G (—oo,+oo) and performing some algebraic manipulations, leads to 




dx [{ux 


-1=0, 4 = 


dxu*^ 


( 6 ) 


Note that for ^ = 1, the integrand of the last term in eq. ® i is a complete differential, which when integrated 
vanishes by the specified asymptotic conditions, hence dl\/dt = 0. For k> 1, to establish conservation laws, 
the goal is to manipulate the integrand of the last term in eq. into a complete differential. Performing 
integration by parts twice, it is readily shown that, for k = 2, the term can be manipulated into a complete 
differential, hence d/ 2 /dt = 0. Meanwhile, for ^ > 2 (and m 0) it can be shown that such a manipulation 
is impossible, hence no further conserved quantities of the form 4 exist. 

More generally, however, certain global conservation laws follow directly from the Hamiltonian struc¬ 
ture of the proposed hierarchy of equations |[25l . For example, the total wave energy 
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is conserved, i.e., dH/dt = 0. Similarly, the total wave mass 
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is conserved, i.e., dM/dt = 0. Fikewise, the total wave momentum 
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is also conserved, i.e., dP/dt = 0. Notice that M = ^4 and P = — ^4- 

From Noether’s theorem (see, e.g., |[26l §4]), these three physical conservation laws arise from un¬ 
derlying symmetries of the governing PDEs, specifically P, PI and M are conserved due fo franslafional 
invariance in space, franslafion invariance in time and fhe abilify fo add arbifrary consfanfs fo the potential 
(p, respectively, as can be easily verified from eq. ([Bl. 

However, besides fhe special case {n,ni) = (1,1) (KdV), we do nof expecf members of fhis hierarchy 
fo be “fully integrable” (i.e., fo possess an infinife number of conservafion laws), yef this remains an open 
question. An important step in this direction would be to establish a Miura transformation 0 between the 
(l,»i) and {2,ni) equations in the proposed hierarchy for m> 1. 
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2.2. Scaling properties 

Above, we have implicitly assumed that we are working with dimensionless quantities, hence no coeffi¬ 
cients involving typical physical parameters appear in the equations. Consider the scaling transformation 

x^xji^ u'^ujW^ ( 10 ) 

where I, z and V are constants, which converts eq. (|5l) from its canonical form to a dimensional form: 

n, + £n% + 5[(n,)X = 0> e = il{rV"), 5 = (11) 

Clearly, given a non-canonical form of one of the equations of the hierarchy (with advective and dispersive 
coefficients e and 5, respectively), we can always choose characteristic scales i, z and V such that the 
rescaling given in eq. (fTOb transforms the PDE (fTT]) i into its canonical given form in eq. ([5]l. 


3. Connection to other nonlinearly dispersive KdV-like equations from the literature 

Clearly, eq. Q reduces to the KdV equation |0 eq. (1)] for {n,m) = (1,1), while it reduces to the 
modified KdV equation |0 eq. (2)] for {n,m) = (2,1). Meanwhile, for {n,m) = (1,3), eq. (I5]l takes the form 
of the unidirectional weakly-nonlinear equation for acoustic waves in an inviscid, non-thermally conducting 
compressible fluid with a material characteristic length coefficient that is a quadratic function of the velocity 
gradient, which was recently derived by Jordan and Saccomandi |[20l . And, for (?i,ni) = (2,3), eq. (I5]l takes 
the form of the unidirectional weakly-nonlinear equation for shear waves in a hyperelastic, incompressible 
solid with material dispersion (of the type suggested by Rubin et al. ll22ll l. which was derived by Destrade 
and Saccomandi ifT^ . The reader is referred to Maugin’s historical review ifTSl for a discussion of various 
other KdV-like equations that arise in the theory of wave propagation in elastic solids. 

Previously, Rosenau and Hyman ifTOll proposed a different hierarchy of nonlinearly dispersive general¬ 
ized KdV equations, namely the K{n,m) equations: 

Ut + (n")jf -|- {u'”)xxx = 0. (12) 

However, the K{n,m) hierarchy is not Hamiltonian, thus lacking a physically desirable property that pro¬ 
vides underlying geometric structure to many nonlinear evolution equations ll27l . A slight modification of 
the K{n,m) hierarchy restores the Hamiltonian structure as shown by Cooper et al. ifTTll . yielding theK*{I,p) 
hierarchy (with certain signs changes and setting a = 1/2 in ifTTl eq. (4)] for consistency of notation): 

u, + u’^'^Ux + uPuxxx + 2puP^'^UxUxx + \p{p-'^)ii’’^^{ufiy’ =0. (13) 

The corresponding Lagrangian density is + /(/ii) — \{(px)^{tPxx)^, which, when compared 

to eq. O, suggests that the present nonlinearly dispersive generalized KdV hierarchy is, in a sense, the 
“most straightforward” generalization of KdV’s Hamiltonian structure. In keeping with the convention in 
the literature, we denote the hierarchy of equations ([51) as K^{n,m). 


4. Traveling wave solutions for the hierarchy 

Introducing the traveling wave ansatz u{x,t) = U{^), where ^ = x — ct for some wave speed c > 0, 
reduces eq. (I5]l to an ordinary differential equation (ODE): 

-cU' + U’V'+[{U')”f = 0, 


(14) 
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where primes denote differentiation with respect to . A first integral of eq. (fTdl) is evident by inspection, 
then rearranging the latter, multiplying by U' and rewriting all terms as complete derivatives immediately 
yields a second integral: 
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where Ci and C 2 are constants of integration. Finally, eq. (fTSl) can be solved implicitly: 
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where we choose the positive root on the left-hand side of eq. (IThl) 1 , hence the ± sign on the right-hand side. 

A complete study of eq. (IThl) for various boundary conditions is beyond the scope of the present work. 
Nevertheless, for (n, m) = (1,1) (KdV), solutions can be written explicitly using hyperbolic or elliptic special 
functions ||4j|5l. Some solutions are also known for m = 3, and an intriguing aspect of the latter is that they 
are compact, i.e., u vanishes outside of some open v-interval. Following |[T^l20]| . let us restrict to localized 
waveforms such that 0 as \x\ -A ± 00 , then Ci = C 2 = 0. It follows that the denominator of the 

integrand on the left-hand side of eq. (fT^ i vanishes atU = { 0 ,[/max}^ where [/max = [(?^ + l)(n-|- 2 )c/ 2 ]'(". 
The Lipschitz condition is not satisfied af fhese poinfs and differenf solutions can be pieced together there, 
allowing for compactification of the waveform. Notice that [/ = 0 is always a solution of eq. (031) for C 2 = 0. 

Now, recognizing the integral on the left-hand side of eq. (fThll i (with Ci = C 2 = 0) as a hypergeometric 
function Il28l §16.15], we obtain 
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where the last integration constant ^0 is fixed by requiring fhaf U (0) = U„ 
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As in MM, if is now sfraighfforward fo piece fogefher a compacf fraveling wave solufion such thaf 
[/ = 0 for G —^ 0 ] U [+^Q, +°°), U is given implicifly by eqs. ([TTIl-lfTS]) for ^ G {—^ 0 , -|-i^o)\{0} and 

C(0) = [/max. 

If can be verified thaf U and U' are confinuous af (§ = {0,±l§o}^ while 


U" = U-^ ([/«,, - [/")-^ ic{n + 1 ) 


t^")(t^max) '"'+'{n + 2)c/{2m), 


(19) 


which suggesf fhere are several possibilifies for U”. The ferm U "'+1 dominates fhe behavior of U" as 
[/ ^> 0+ ( 1 ^ —^ ±i§o) wifh U” -A 0 for m < 3, U” -A y/c/6 for m = 3 (see also ||20l eq. (4.18)]) and |[/"| —> 00 

for m > 3. Similarly, fhe ferm ([/^ax “ [/")^^ dominates the behavior of U” as U ^ U max 0) and, 
since the compact solutions derived hold only for m > 1 and n > 0, | U" I —>■ 00 as [/ —)■ U max 0 )- 

The presence of an infinite second derivative at = 0 means that the compact waveforms constructed are 
peaked like the (non-compact) peakon solutions of the Camassa-Holm equation 0, hence we refer to the 
family of solutions given by eqs. (fT7]) - (fT8]) by the portmanteau peakompactons. The possibility of an infinite 
second derivative as —)■ { 0 , ± 1 ^ 0 } implies that we have constructed pseudo-classical solutions of eq. ([S]), 
i.e., piecewise solutions that satisfy eq. ([5]) in a classical sense away from ^ = {0,±(^o}^ and, moreover, 
lead to well defined limifs for fhe nonlinearly dispersive terms m{m — l)([/')'"~* {U'')^ and m{U')’^^^U"' as 
^ -A {0,=h(§o} (see, e.g., fhe illuminafing discussion in 1291 |30l). 

The peakompacfon solufion (fTTll - dTS]) is illusfrafed in fig.[T]for various values of {n, m). Nofe fhaf Umax is 
independenf of m, hence fhe form of fhe nonlinear dispersion has no effecf on fhe peak heighf. However, both 
n and m affect the length of the finite support since (§0 is a function of both. Similarly, c is featured in both 
the expression for Umax and (§o^ thus the wave speed affects both the width and height of the peakompactons. 
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Fig. 1. Illustration of the K*(n,m) hierarchy’s peakompacton solutions given by eqs. (nJ-dlU for c = 0.75 and (a) n = 1,2,3 
(solid, dashed, dotted) and m = 3, (b) n = 1 and ra = 3,4,5 (solid, dashed, dotted). 


5, Conclusion 

In the present work, we proposed a hierarchy of nonlinearly dispersive generalized Korteweg-de Vries 
(KdV) evolution equations based on a modification of the Lagrangian density whose induced action func¬ 
tional is extremized by the KdV equation. Conservation laws, scaling properties and peaked, compact trav¬ 
eling wave solutions for this class of PDEs were discussed. It is important to emphasize that the proposed 
K^{n,m) hierarchy of nonlinearly dispersive KdV-like PDEs possesses compact traveling wave solutions, 
like the K{n,m) hierarchy. However, unlike the K{n,m) hierarchy, the K^{n,m) equations also possess 
Hamiltonian structure. Indeed, the K*{l,p) hierarchy restored the Hamiltonian structure to the K{n,m) 
equations, however, to the best of our knowledge, none of the equations from the K*{l,p) hierarchy have 
been derived from physical principles, i.e., these equations remain phenomenological. On the other hand, 
at least two of the equations from the proposed K^{n,m) hierarchy (with m > 1) have been independently 
derived to describe wave propagation in generalized continua. 

In the present work, our goal was to present the K*{n,m) class of nonlinear evolution equations and 
some of their elementary properties. The mathematical analysis of these PDEs remains to be done. A 
number of questions were noted in the text above: e.g., Miura transformations between members of the class, 
integrability of the hierarchy and/or individual members of it, existence of other (exotic) exact solutions, 
collision properties of the known solitary waves, etc. However, of specific interest would be to understand 
whether eq. dSll with m even (e.g., 2) has a physical interpretation, since (so far) only eq. dSll with m odd 
(specifically, 1 and 3) has been derived in physical contexts. 


Acknowledgments 

I.C.C. was supported by the LANL/LDRD Program through a Eeynman Distinguished Eellowship. Los 
Alamos National Laboratory (LANE) is operated by Los Alamos National Security, L.L.C. for the Na¬ 
tional Nuclear Security Administration of the US Department of Energy under Contract No. DE-AC52- 
06NA25396. Prof. Andrus Salupere and the local organizing and international scientific commitfees of fhe 
2014 lUTAM Symposium on Complexity of Nonlinear Waves in Tallinn, Estonia are acknowledged for their 
work in making the meeting a great success. The symposium’s hospitable intellectual environment made 
the present work possible. Einally, I.C.C. would also like to thank E. Cooper, J.M. Hyman, PM. Jordan, T. 
Kress, A. Oron and A. Saxena for helpful discussions on the topic of the present work. 


REFERENCES 

1. J. Boussinesq. Essai sur la theorie des eaux courantes. Memoires presenles par divers savants a VAcademic des Sciences de 
rinstitutde France, XXIII:l-680, 1877. 








6 


Proceedings of the Estonian Academy of Sciences, 


2. D. J. Korteweg and G. de Vries. On the change of form of long waves advancing in a rectangular canal, and on a new type of 

long stationary waves. Phil. Mag. (Ser. 5), 39:422-443, 1895. 

3. N. J. Zahusky and M. D. Kruskal. Interaction of solitons in a collisionless plasma and the recurrence of initial states. Phys. 

Rev. Lett., 15:240-243, 1965. 

4. R. M. Miura. The Korteweg-de Vries equation: A survey of results. SIAM Rev, 18:412^59, 1976. 

5. A. Jeffrey and T. Kakutani. Weak nonlinear dispersive waves: A discussion centered around the Korteweg-de Vries equation. 

SIAM Rev, 14:582-643, 1972. 

6. R. M. Miura. Korteweg-de Vries equation and generalizations. I. A remarkable explicit nonlinear transformation. J. Math. 

R/iyi., 9:1202-1204, 1968. 

7. R. S. Johnson. On the inverse scattering transform, the cylindrical Korteweg-De Vries equation and similarity solutions. Phys. 

Left. A, 72:197-199, 1979. 

8. T. R. Marchant and N. F. Smyth, The extended Korteweg-de Vries equation and the resonant flow of a fluid over topography. 

J. Fluid Mech., 221:263-287. 

9. R. Camassa and D. D. Holm. An integrable shallow water equation with peaked solitons. Phys. Rev Lett., 71:1661-1664, 

1993. 

10. R Rosenau and J. M. Hyman. Compactons: Solitons with finite wavelength. Phys. Rev Lett., 70:564-567, 1993. 

11. F. Cooper, H. Shepard and P. Sodano. Solitary waves in a class of generalized Korteweg-de Vries equations. Phys. Rev E, 

48:4027^032, 1993. 

12. C. I. Christov and M. G. Velarde. Dissipative solitons. Physica D, 86:323-347, 1995. 

13. A. S. Fokas. On a class of physically important integrable equations. Physica D, 87:145-150, 1995. 

14. A. Salupere, J. Engelbrecht and G. A. Maugin. Solitonic structures in KdV-based higher-order systems. Wave Motion, 34:51- 

61, 2001. 

15. G. A. Maugin. Solitons in elastic solids (1938-2010). Mech. Res. Commun., 38:341-349, 2011. 

16. J. Engelbrecht and E. Pastrone. Waves in microstructured solids with strong nonlinearities in microscale. Proc. Estonian Acad. 

Sci. Phys. Math., 52:12-20, 2003. 

17. M. Randriitit, A. Salupere and J. Engelbrecht. On modelling wave motion in microstructured solids. Proc. Estonian Acad. Sci., 

58:241-246, 2009. 

18. A. Salupere, M. Lints, J. Engelbrecht. On solitons in media modelled by the hierarchical KdV equation. Arch. Appl. Mech., 

9-11:1583-1593, 2014. 

19. M. Destrade and G. Saccomandi. Solitary and compact-like shear waves in the bulk of solids. Phys. Rev E, 73:065604, 2006. 

20. P. M. Jordan and G. Saccomandi. Compact acoustic travelling waves in a class of fluids with nonlinear material dispersion. 

Proc. R. Soc. A, 468:3441-3457, 2012. 

21. P. Rosenau and A. Oron. On compactons induced by a non-convex convection. Commun. Nonlinear Sci. Numer. Simulat., 

19:1329-1337, 2014. 

22. M. B. Rubin, P. Rosenau and O. Gottlieb. Continuum model of dispersion caused by an inherent material characteristic length. 

J. Appl. Phys., 77:4054-4063, 1995. 

23. C. S. Gardner. Korteweg-de Vries equation and generalizations. IV. The Korteweg-de Vries equation as a Hamiltonian system. 

J.Math. Phys., 12:1548-1551, 1971. 

24. I. M. Gelfand and S. V. Fomin. Calculus of Variations. Dover Publications, Mineola, NY, 2000. 

25. G. A. Maugin and C. I. Christov. Nonlinear Duality Between Elastic Waves and Quasi-particles. In C. I. Christov and A. Guran, 

Selected Topics in Nonlinear Wave Mechanics, pages 117-160. Birkhauser, Boston, 2002. 

26. M. D. Kruskal and N. J. Zahusky. Exact invariants for a class of nonlinear wave equations. J. Math. Phys., 7:1256-1267, 1966. 

27. D. D. Holm, T. Schmah and C. Stoica. Geometric Mechanics and Symmetry: From Finite to Infinite Dimensions. Oxford 

University Press, New York, 2009. 

28. A. B. Olde Daalhuis. Hypergeometric Eunction. In F. W. J. Olver, D. W. Lozier, R. F. Boisvert and C. W. Clark, editors, NIST 

Digital Library of Mathematical Functions, Chapter 15, http: //dlmf . nist . gov/ , release 1.0.9, 2014. 

29. Y. A. Li and P. J. Olver. Convergence of solitary-wave solutions in a perturbed bi-Hamiltonian dynamical system 1. Com¬ 

pactions andpeakons. Discr. Contin. Dyn. Syst. A, 3:419^32, 1997. 

30. Y. A. Li and P. J. Olver. Convergence of solitary-wave solutions in a perturbed bi-hamiltonian dynamical system 11. Complex 

analytic behavior and convergence to non-analytic solutions Discr. Contin. Dyn. Syst. A, 4:159-191, 1998. 


